Tuning a Resonance in the Fock Space: Optimization of Phonon Emission in a 

Resonant Tunneling Device 
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Phonon-assisted tunneling in a double barrier resonant tunneling device can be seen as a resonance in 
the electron-phonon Fock space which is tuned by the applied voltage. We show that the geometrical 
parameters can induce a symmetry condition in this space that can strongly enhance the emission 
of longitudinal optical phonons. For devices with thin emitter barriers this is achieved by a wider 
collector barrier. 
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Progress in mesoscopic semiconductor devices and 
molecular electronics is driven by the need of minia- 
turization and the wealth of new physics provided by 
coherent quantum phenomena. A fundamental idea be- 
hind these advances was Landauer's view that conduc- 
tance is transmittance H,^. Hence, the typical conduc- 
tance peaks and valleys, observed when some control pa- 
rameter is changed, are seen as fringes in an interfer- 
ometer. However, the many-body electron-electron (e- 
e) and electron-phonon (e-ph) interactions restrict the 
use of this picture. The e-e effects received much at- 
tention in different contexts ||]. While interest on e-ph 
interaction remained mainly focused on double barrier 
Resonant Tunneling Devices (RTD) where phonon- 
assisted tunneling shows up as a satellite peak in a valley 
of the current-voltage (l-V) curve, recent observation of 
electro-mechanical effects in molecular electronics Q re- 
quires a reconsideration of the e-ph problem. Theory 
evolved from a many-body Green's function in a simpli- 
fied model for the polaronic states |^] to quantum and 
classical rate equations approach The latter uses an 
incoherent description of the e-ph interaction by adopt- 
ing an imaginary self-energy correction to the electronic 
states 

In this article, we analyze a quantum coherent solu- 
tion of transport with e-ph interaction. We resort to 
a mapping of the many-body problem into a one-body 
scattering system where each phonon mode adds a new 
dimension to the electronic variable |ll|Jl2t] . Transmis- 
sion of electrons between incoming and outgoing channels 
with different number of phonons are then used in a Lan- 
dauer's picture where the only incoherent processes occur 
inside the electrodes. This allows to develop the concept 
of resonance in the e-ph Fock space and the identification 
of the control parameters that optimize the coherent pro- 
cesses leading to a maximized phonon emission. It also 
gives a clue as to how "decoherence" arise within an exact 
many-body description. As an application, we consider a 



RTD phonon emitter where the relevant parameters are 
best known. There, the first polaronic excitation serves 
as an "intermediate" state for the phonon emission. An 
electron with kinetic energy e < ep and potential en- 
ergy eV in the emitter decays through tunneling into an 
electron with energy e -I- e\/—hujQ in the collector plus a 
longitudinal optical (LO) phonon. The tuning parameter 
is the applied voltage while the optimization of phonon 
emission requires the tailoring of the tunneling rates. 
Let us consider a minimal Hamiltonian: 
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(1) 



where c^and c,- are electron creation and annihilation 
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operators at site j on a 1-d chain with lattice constant 
a and hopping parameters Vjj+i — V. Tunneling rates 
are fixed by Vb,i = Vr and V-i^ = Vl (14(r) < V). 
The site energies are Ej = 2V for j < and 2V — eV 
for j > 0. Eq ~ -B(o) — a eV is the well's ground state 
(including the charging effect) shifted by the electric field. 
For barrier widths Ll and and well size Lw a linear 
approximation for the potential profile gives a = {Li^ + 
Lw/2)/(LL-f -Lw + -^r). We consider a single LO-phonon 
mode and an interaction (Vg) limited to the well. 6+ and 
b are the creation and annihilation operators for phonons. 
We restrict the Fock space to that expanded by \j, n) = 
(6+)" /\/nI |0) , which maps to the 2-dimensional one- 



body problem shown in Fig. 



The number n of phonons 



is the vertical dimension 12| . The horizontal dangling 
chains can be eliminated through a decimation procedure 
]l0|,p3[ leading to an effective Hamiltonian: 



Ue-ph - ^{[^^0 + nhuja + E„(e)] |0, n) (0, n\ 



(2) 



-^/W+\Vg (|0, n + 1) (0, n\ + |0, n) (0, n + 1|)}, 

The electron hopping into the electrodes is taken into ac- 
count by the £-dependence of the retarded self-energy 
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FIG. 1. Simple model: dots are states in the Fock space, 
lines are interactions. The efltective Hamiltonian including 
two entangled polaronic states is represented at the bottom. 



corrections !]„ — + ^S]„ . Specifically, = 



with: 



X i:{e-nhuJo), = |-^| x Y. {e - nhcoo + eV) , 



e — e' 



I](e) = A(e)-ir(£); A (e) = i j 

r (e) =\lv^- {e/2 + Vf x e{e) x e{AV - e). (3) 

While the imaginary part F = hv^ja is proportional to 
the group velocity in the electrodes, the actual es- 
cape rates rL(R)/?i are barrier controlled. For width 

I 1 2 

Ll(r) and attenuation length ^ , rL(R)/F = |Vl(r)/V| 

~ exp[-LL(R)/C]- 

The retarded Green function connecting states i and 



G^,, (e) = (0, n\ {eX - He-ph (e))"' |0, i) 



(4) 



has poles at the exact eigenenergies. If = 0, these 

are the polaronic energies Eq — + nhujQ. The trans- 
mission coefficient Tnj from the i-th incoming channel at 
left electrode to the n-th channel at right is |lO| : 

r„,,(e)^2ImpE„(£)] \G^^^{ef 2Im[^Y,ie)]. (5) 

If the Fermi energy ep ^V, Eq.(|^) becomes: 

E(e)«-ir(e = ei.) x0(e), (6) 

and the 0-function may cancel some T's. 

To obtain the elastic transmittance when g = 
{Vg/hivoy <C 1 and (e_F,rL -t- Fr) < hu>o, we need: 
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0,0 



e-En 
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ipL + Tr] 
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e-[Eo + Tiuo] + ipL + Tr] ' 



(7) 



evaluated with the first two polaronic states. Here, 
Fl = 5Fl and Eq ^ Eq - -1^. The first term con- 
tains the main resonance associated to the build up of 
the polaronic ground state. The second term contains 
a virtual exploration into the first polaronic excitation. 
It is noteworthy that when F = 0, this Green function 
would cancel out at an intermediate energy giving rise 
to an antiresonance p^ , p^ . This concept extends the 
spectroscopic Fano- resonances [|5| to the problem of con- 
ductance For g <C 1, this effect is less important 
and in the whole energy range. 



T, 



4rLFR 



0,0 



[e-So] +[rL + rR]2 



0{g) 



(8) 



describes the main resonant elastic peak at e = Eq 

The inelastic transmittance 
from: 



Ti Q, can be evaluated 



<-^i,o 



'0 



En 



i[rL + Tr] 

hujo 



+ 



e-[Eo + huo] + i[FL + Fr] 



(9) 



When e + eV > hujo escapes are enabled and its poles 
involve the processes represented in the inset of Fig. ||: 
a) The first term gives an inelastic transmittance at the 
main peak, b) The second term provides a satellite peak 
at e = Eq + tkoo, associated to a polaronic excitation fol- 
lowed by its decay into an escaping electron and a phonon 
left behind. Around this satellite peak: 
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[e - {Eo + hcjo)] +[rL + FR] = 



(10) 



showing that phonon emission is a resonance in the Fock- 
space (see bottom of Fig. |l|). A maximal probability 
(Ti Q = 1) requires equal rates of formation and decay 
|]l^: Fl = Fr, which in our RTD imphes: 



2Cln[ 



V„ ^ 



(11) 



Hence, thin barriers with this generalized symmetry con- 
dition have T]^ Q ~ 1 over a broad energy range. 

The application of the Keldysh formalism to our 
Fock-space gives an electrical current l^ot expressed as a 
balance equation Q| in terms of the transmittances of 
Eq. (||) and the electrochemical potentials. The ex- 
perimental condition of high bias and low temperature 
(eV > ep ^ fcsT), rules out right-to-left flow, while 
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FIG. 2. Current density as a function of the applied voltage 
for a symmetrical (thin line) and optimized (thick line) struc- 
tures with Ll = 7a (19.7 A). The dotted line indicates the 
background current in the region of the satellite peak for the 
symmetrical structure. The inelastic processes contributing 
to the peaks a) and b) are represented in the inset. 

phonon production should maximize the emitted power 
according to Eq. (pl|). 

Let us compare these basic predictions with the numer- 
ical results of a description involving geometry, voltage 
and energy dependences of a typical RTD. A discrete 3-d 
model is defined in terms of the effective mass m* with 
V = / {2m* a?) . The potential profile for the diagonal 
energies Ej is shown in the inset of Fig. Ni,, TVr and 
iVw are the number of sites in the left and right barri- 
ers and the well, the associated widths are Li = Nia. 
For translational symmetry along the interface, we con- 
sider a single phonon mode per transversal (parallel to 
the interface) state, with frequency wq and localized in 
the structure region. While conservation of transverse 
electron's momentum might not be fully realistic it 
constitutes a first approximation yielding results consis- 
tent with the main experimental features. The current 
components are obtained from ( p^ ) by integration over 
the transversal modes. The parameters in our calcula- 
tions are chosen to simulate the case of a GaAs-AlGaAs 
structure. The effective mass m* is 0.067 me, the LO 
phonon frequency Tluiq = 36 meV, a = 2.825 A, and the 
hopping parameter V = 7.125 eV. A typical e-ph interac- 
tion strength of g ~ 0.1 is obtained with ~ 10 meV. 
For a well of 56. 5A, barrier heights of 300 meV and £f = 
10 meV, the inclusion of rt < 3 warrants good numerical 
convergence. 

For wide left harriers (of about 25a ~ 70 A or more), we 
found that the maximum value of Pin varies slowly with 
the width of the right barrier. Hence, consistently with 
our discussion of the 1-d model, there is no substantial 
gain in Pjn by choosing an asymmetric structure. Con- 
sequently, a high phonon emission rate should be sought 
for thin barriers. 
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would require a 54A collector's barrier. We expect that 
our results could stimulate the search for excited phonon 
modes (e.g. with Raman spectroscopy), in operational 
RTD's as a function of the applied voltage in the various 
configurations. While for simplicity we have restricted 
our analysis to a model RTD, our analysis applies to 
other problems involving electronic resonant tunnel- 
ing in the presence of an interaction with an elementary 
excitation. 

We acknowledge financial support from CONICET, 
SeCyT-UNC, ANPCyT and Andes- Vitae-Antorchas. 
HMP and LEFFT are afhhated with CONICET. 



FIG. 3. Power emitted as LO phonons as a function of the 
applied voltage for Ll ~ 7a (19.7 A) and different values of 
Lr. The efficiency as a function of the right barrier width is 
shown in the inset. 



A thin left barrier of Ll = 7a (19. TA), gives a tunneling 
probability Ti^iep) = Pl/F - 0.03. Figure| shows the I- 
V curves for symmetric and asymmetric RTDs. In Fig. ^ 
we show Pin— V for various right barrier widths Lr. The 
peaks are shifted to higher voltages as Lr is increased, 
because the resonant energies are lowered approximately 
by aeV. We can also see that the peak value of Pin as a 
function of the right barrier width exhibits a maximum. 
The l-V curve for the optimal configuration is shown in 
Fig. H (heavy line). The inset of Fig. || shows the de- 
pendence of ry, evaluated at the optimal voltage, on the 
right barrier width for various left barrier widths where 
Pl < huQ. In agreement with our theoretical analysis, rj 
keeps the same magnitude for all the shown geometrical 
configurations. The main result of Fig. ^ is the con- 
firmation that, for a given Ll satisfying Pl < hujo a-nd 
gPL + Pr > e_F, the phonon emission rate is enhanced 
by a factor 2.5 by choosing a wider right barrier as pre- 
scribed by Eq. (|l^) . This may explain the unusually large 
satellite peaks of asymmetric structures [|l^ . 

An RTD optimized for phonon emission might have 
many applications. In fact, in AlGaAs-GaAs RTD these 
primary LO phonons have a short life-time ||20|| and decay 
into a pair of LO and transverse acoustic (TA) phonons. 
This phenomenon inspired the proposal | [Tl| | for the gener- 
ation of a coherent TA-phonon beam in an RTD (called 
a SASER) |^l|]. That device required an energy differ- 
ence between the first two electronic states in the well 
El 



Eq — Tiloq 1 11 2l|]. In contrast, the present pro- 
posal does not require such an accurate device geometry. 
Instead, operation in the phonon emission mode only re- 
quires the tuning of the many-body resonance with the 
external voltage. Geometry just improves its yield by 
imposing a generalized symmetry condition in the Fock- 
space. For a typical AlGaAs emitter barrier of 20A this 
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